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Motivation 


Partition function : 
Z = tr[exp -6H (p, ĉ)] 


Path integral representation 


B 
Z= | DeDpexp | dt (ipt — H(p,x)). 
0 


Infinitely multiple integration! 
e@ Mean field theory + perturbations 


ə Monte Carlo sampling 
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QCD phase diagram(?) 


Temperature 7 


Sj Color superconductors) 
CFL-K®, Crystalline CSC 


Meson supercurrent Baryon Chemical Potential us 
Gluonic phase, Mixed phase 


Nuclear Superfluid 


(Fukushima, Hatsuda, 2010) 


B 
Z(6,u)= foa det(y;D;(u, A)) exp er d’xtrF? 


DAG 


Lefschetz-thimble path integral 


Motivation 


Path integral with complex weights appear in many important 
physics: 


e Finite-density lattice QCD, 
spin-imbalanced nonrelativistic fermions 


ə Gauge theories with topological 6 terms 


ə Real-time quantum mechanics 


Oscillatory nature hides many important properties of partition 
functions. 
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Example: Airy integral 
Let's consider a one-dimensional oscillatory integration: 


a= | i T as 
p 3 - 


RHS is well defined only if Ima = 0, though Ai(z) is holomorphic. 


Horo nap 


exp( osx ested) 


AAACN 


Figure : Real parts of integrands fora = 1 (x10) & a = 0.5i 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Feb 20, 2015 @ Keio U. 6 / 38 


Contents 


ə How can we circumvent such oscillatory integrations? 
=> Lefschetz-thimble integrations 
[Witten, arXiv:1001.2933, 1009.6032] 


ə Applications of this new technique for path integrals 


> Phase transition of the matrix model 

[YT, Phys. Rev. D91, 036002 (2015); Kanazawa, YT, arXiv:1412.2802] 
> Quantum tunneling in the real-time path integral 

[YT, Koike, Ann. Physics 351 (2014) 250] 
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Lefschetz-thimble integrations 


Introduction to Lefschetz-thimble integrations 
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Again Airy integral 


Airy integral: 


3 
Ai(a) = f T epi (5 + ar) À 


The integrand is holomorphic w.r.t x 
=> The contour can be deformed continuously without changing the 
value of the integration! 
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Lefschetz-thimble integrations 


Steepest descent contours 
What is the most appropriate contour for our purpose? 


Re[iS(a,a)] should be made as small as possible. 
<> The contour should be perpendicular to Re[iS(, a)] = cosnt. 


Steepest descent ones J must satisfy Im[iS(x,a)] = const. 
because of the holomorphy. 


m 


5 i F i, | 
m. Contour plots for ReliS(x, a)] with a= = exp u 
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Lefschetz-thimble path integral 


Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (o = 1,2) for the Airy 


integral: 
d 3 
Ai(a) = > No [ Pi (3 taz) 


Ng: intersection number of the steepest ascent contour Ke and R. 


il 
| 
} 
+ 
1 
| 
i 
| 


STEHE 

| 
| | 
3 


Figure : Lefschetz thimbles J and duals K (a = exp(0.1i), exp(7i)) 


3 
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Airy integral & Airy equation 


Let us consider the “equation of motion”: 


dz d eilz*/3+az) _ Q, 
Qn dz 


This is nothing but the Airy equation: 


Two possible integration contours = Two linearly independent 
solutions of eom. 
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Lefschetz-thimble integrations 


Generalization to multiple integrals 


Model integral: 


z= | dzı -- - d£n exp Z(x;). 


What properties are required for Lefschetz thimbles J? 
@ J should be a n-dimensional object in C”. 
© Im[Z] should be constant on 7. 
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Lefschetz-thimble integrations 


Short note on technical aspects 
Complexified variables (a = 1,..., n): Za = Za + ipa. 
Regard x, as coordinates and p, as momenta, so that Poisson 
bracket is given by 


N [seat - so 


a=1,2 


O1 Opa OLa Opal 
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Lefschetz-thimble integrations 


Short note on technical aspects 
Complexified variables (a = 1,..., n): Za = £a + ipa. 
Regard x, as coordinates and p, as momenta, so that Poisson 
bracket is given by 


u Of Og ög Of 
gh = 2 Ee Ops: On, an 


Hamilton equation with the Hamiltonian H = Im[Z(z,)]: 


df (x, p) dza OT 
ig eee (= dt --(#)) 


This is Morse’s flow equation: “Re[Z] < 0. 


= We can find n good directions for J around saddle points! 
[Witten, 2010] 
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Multiple integral: Lefschetz-thimble method 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles Js: 


d?r eiS(2) = ne | d"z ee, 
I 2 Jy 


Js are called Lefschetz thimbles, and Im[iS] is constant on it. 


Ne. intersection numbers of duals X, and R”. 
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Searle E EEC Harmonic oscillator 


Example of QM: Harmonic oscillator 


The action functional 


Saddle-point condition (=Euler-Lagrange eq.): 


d?z.ı(t 
zal ) = — zalt). 


The classical solution: 
Te — T; cos T 


zalt) = ao ae sin t + x; cost. 


Lefschetz-thimble path integral 
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bent caeli 
Example of QM: Flow equation 


Flow equation 
o o? — 


Boundary conditions: Az(t; —co) = 0, and Az(0; u) = Az(T, u) = 0. 
The set of solutions Az are spanned by 


eit/4 exp (ey = 1) u sin", (nr >T), 
ale pdas. men 


Azn(t; u) = 
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Penit eee 
Example of QM: Lefschetz thimble 


Lefschetz thimble J (~ R”): 


. net = wet 
—im/4 Bee im/4 MR 
za(t) + € X usin T +e ` arsin = | ae € R 
é=1 f=v4+1 
v: the maximal non-negative integer smaller than T/r. 
Integration measure on J becomes 
V . (oe) . (oe) 
/ Dz = nf I] e idan I] ei dan = eeN |T] Vidap. 
I n=1 m=v+1 f=1 
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binet sa. 
Example of QM: Path integral 


Feynman kernel for the harmonic oscillator: 


| 
© 


TZ 
E Zei] _ N vin E 
Tze) _ ‚TU 1 1 
xp (A E / Var jl T/r] 
= Alea oe us 
u 2rihl aA 2 


| 
© 


Appearance of the Maslov-Morse index v is quite manifest. 
[YT, Koike, Ann. Physics 351 (2014) 250] 
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Haven oral es 
Proposed Applications of Lefschetz-thimble methods 


ə Analytic properties of partition function for CS theories 
[Witten] 


ə Real-time description of tunneling 
[YT, Koike; Cherman, Unsal] 


oe Sign problem of the lattice Monte Carlo simulation 


[M. Cristoforetti, et al., Kikukawa, et al.] 


e@ Geometrization of resurgence transseries 


[Basar, Dunne, Unsal; Cherman, Dorigoni, Unsal] 
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Gross-Neveu-like model 


Phase transition of a matrix model 
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Gross-Neveu-like model 


0-dim. Gross-Neveu-like model 


The partition function of our model study is the following: 


Zu(G,m) = | apaw exp [Sua + mya t ur ( ya) | l 


The Hubbard-Stratonovich transformation gives 


N 
Zn(G,m)=4/ “a faces, 


— log[p* + (o + m)’]. 


with 


For simplicity, we put m = 0 in the following. 
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Gross-Neveu-like model 


Properties of S(c) 
S has three saddle points: 


_ OS(z) _ 22 2z 


0 u 
Oz G P+2 


= z=0 +/G-P. 


Figures for G = 0.7eF"" and p? = 1 [Kanazawa, YT, arXiv:1412.2802]: 
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Figures for G = 1.1eF®"! and p? = 1 [Kanazawa, YT, arXiv:1412.2802]: 
8 P 
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From these figures, we learn that, for real G, 


ə z = is the unique critical point contributing to Z if G < p°. 
ə All three critical points contribute to Z if G > p°. 
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Gross-Neveu-like model 


Stokes phenomenon 


The difference of the way of contribution can be described by 
Stokes phenomenon. 

< At some special values of coupling, several critical points are 
connected by the flow. [Witten, arXiv:1001.2933, 1009.6032] 


Figures of G-plane for ImS(0) = ImS (z4) [Kanazawa, YT, arXiv:1412.2802]: 


-1.5 


“10-05 00 05 10 15 20 


RG 
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Gross-Neveu-like model 


Dominance of contribution 


The Stokes phenomenon tells us the number of Lefschetz thimbles 
contributing to Zy. 


However, it does not tell which thimbles give main contribution. 


In order to obtain (a) # 0 in the large-N limit, z+ should dominate 
z=0. 


> ReS(z;) < ReS(0) 
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Connection with Lee-Yang zero 


Blue line: ImS(z+) = ImS(0). 
Green line: ReS(z+) = ReS(0). 
Red points: Lee-Yang zeros at N = 40. [Kanazawa, YT, arXiv:1412.2802] 
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Gross-Neveu-like model 


Conclusions for studies with GN-like models 


@ Decomposition of the integration path in terms of Lefschetz 
thimbles is useful to visualize different phases. 


@ The possible link between thimble decomposition and Lee-Yang 
zeros is indicated. 
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Resergence theory 


Resurgence theory of the perturbation series 
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Perturbation theory 


Let us consider the integral 


Z(d) = [aves l-4 (= = >) | 


Gaussian approximation around x = +1 gives 


ZA) ~ ANY, anA” 


n=0 


~~ TT 2m — 1)? — 1/4) 


(an is growing factorially in terms of n. ) 
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Resergence theory 


Non-convergent series 


M-test says the convergence radius X. is 


An n+1 
As < = — 
u An+1 (2n + 1)? = 1/4 


0. 


Perturbation series = Asymptotic series (Non-convergent!) 


Can we make sense of it? 
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Resergence theory 


Borel summation 


Define 


This is convergent at least for small t. Since 


f dit?e > = nlart!, 
0 


we might be able to consider the following series 


el/4r oo 


= = dtBP(t)e" 
TE 


BA) 


However, BP(t) can have a pole on R4. 
= This integral becomes ambiguous. 
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Resergence theory 


Interpreting Borel ambiguity via Lefschetz thimbles 


Ambiguous imaginary part of the Borel sum = Jo is missed!. 
(Resurgence trans-series theory: Ecalle, Voros, etc. ) 
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Lesson & application to QFT 


Lesson: 


e Large order perturbation theory can become ambiguous even 
after resummation of the series. 


ə Such ambiguity provides information on contributions of 
non-trivial classical solutions. 


Application to QM & QFT 
ə This provides an exact semiclassical treatment of QM & QFT. 
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Example of nontrivial solutions: Gauge theory 


Diakonov et. al., Phys. Rev. D70, 036003 (2004) 


(KvBLL caloron) 
These are related to dynamics of QCD! (Ünsal, Schaefer, Poppitz, Dünne, 
Yuya Tanizaki (University of Tokyo, RIKEN) 
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Resergence theory 


Example of nontrivial solution: Double-well QM 


Imaginary-time instantons as exact solutions 
(YT, Koike, Ann. Phys. 351 (2014) 250): 


(n,m)=(1,0) (n,m)=(3,0) 
Sp Lol p= 
Y 4 
05 al / | 
2 4 8 io 12 \ 7 
| | | 
-0.5 | -0.5 / / 
/ pei / A — Rea 
en ae zwä aija — =, ER 


4 
Seisssiesl = 3 (n + m) 


ı - = = 
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Resergence theory 


Real-time instanton solutions 


Example of nontrivial solution: Double-well QM 


(vr, Koike, Ann. Phys. 351 (2014) 250; Cherman, Unsal, arXiv:1408.0012): 
(n,m) = (31,30) 


(n,m) = (52,50) 
AAA N 


w 
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Summary 


ə Lefschetz-thimble integral is a useful tool to treat multiple 
integrals. 


ə Semiclassical expansion of QM can be reformulated in a neat 
way with this technique. 


@ Perturbation theory in a strong coupling regime (7?) 


ə Sign problem of Monte Carlo study (7?) 
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